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Abstract 

We perform a general analysis of the contributions to the muon anomalous magnetic 
moment in models with a superlight gravitino. We discuss the interpretation and 
the phenomenological implications of the results. We find that present constraints 
on the model parameters are comparable and complementary to the ones coming 
from collider searches and from perturbative unitarity. However, the Brookhaven 
E821 experiment will probe large unexplored regions of parameter space. 
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1. Introduction 



The anomalous magnetic moment of the muon is both one of the best measured quanti- 
ties in particle physics and one of the most accurately known predictions of the Standard 
Model (SM) of strong and electroweak interactions. In terms of = ((7^ — 2)/2, the 
present experimental and theoretical situation can be summarized [|I| as follows: 



a"-^ = (116592350 ± 730) x 10"" 
> 



> ^ ^ ' (I) 

= (116591596 ± 67) x 10"". ^ ^ 



As a consequence, powerful constraints on extensions of the SM can be extracted from the 
comparison of eq. (|l]) with possible non-standard contributions, 6a^. The 95% confidence 
level bound corresponding to the present data isQ 

-0.7 X 10"^ < < 2.2 X 10"^ (2) 

and the E821 experiment at Brookhaven is expected to reduce further the experimental 
error by roughly a factor of 20. 

The most popular and theoretically motivated extensions of the SM are those incor- 
porating low-energy supersymmetry. The study of in supersymmetric models has a 
long history. Already in the early days of supersymmetry it was realized that = 
in the limit of unbroken supersymmetry, both in the global and in the local 0] case. In 
the case of broken supersymmetry, the first estimates of 6a^ from loop diagrams involving 
supersymmetric particles were given in |]5|. By now, the full one-loop contribution to 
in the Minimal Supersymmetric extension of the Standard Model (MSSM) is available, for 
an arbitrary supersymmetric particle spectrum P, 0: the relevant one-loop diagrams in- 
volve smuon-smuon-neutralino and chargino-chargino-sneutrino loops. Thanks to the fast 
decoupling properties of the soft supersymmetry-breaking mass terms, these contributions 
are comfortably within the bounds of eq. (|^) for typical values of the MSSM parameters 
allowed by the direct searches for supersymmetric particles. The only situation in which 
Sttfj, can saturate the bounds of eq. (0) is when the masses of supersymmetric particles 
are close to their present lower bounds and tan/3 = V2/V1 is very large (so that the muon 
Yukawa coupling is considerably enhanced). 

In supersymmetric models with a light gravitino, the effective low-energy theory does 
not contain only the MSSM states, but also the gravitino and its superpartners, and we 
must compute the additional contributions coming from the one-loop diagrams where 
these particles (as well as other MSSM particles) are exchanged on the internal lines. Most 
of the existing calculations of these effects were performed in the framework of supergravity, 
and the discussion [|, ^, [10], |Tl]] was mainly focused (with some controversies) on the 



finiteness of the broken supergravity contributions. The most complete study available so 



-'^More conservative estimates of the present theoretical error associated with the hadronic contribution 
can enlarge the total theoretical error by up to a factor of two. Given the present experimental precision, 
this would not affect significantly the bounds of eq. (0), and improvements are expected in the future. 
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far is the one of ref. |Ty], which considered a general spontaneously broken supergravity, 
and computed the one-loop diagrams involving virtual particles from the matter, gauge 
and gravitational supermultiplets. However, in the case of a light gravitino (the only 
phenomenologically relevant one for this type of study) we can work directly in the globally 
super symmetric limit, keeping only the goldstino and its spin-0 superpartners (sgoldstinos) 
as the relevant degrees of freedom from the supersymmetry-breaking sector [|12|. Indeed, 
this method has already been partially applied in , to compute the contribution to of 
the goldstino-smuon-smuon loops. In this paper, we shall adopt such a method to perform 
a general analysis of the contributions to 5a^ from the supersymmetry-breaking sector of 
the theory. 

The paper is organized as follows. In sect. 2, we introduce the model and derive the 
relevant interaction terms. In sect. 3, we present the general expressions for our results. In 
sect. 4, we discuss their interpretation, with emphasis on possible divergent contributions. 
In sect. 5, we specialize our formulae to some phenomenologically interesting limits and 
discuss their implications. Whenever possible, we compare with the previous literature, 
in particular with the supergravity computation of ref. |1T0|. 



2. The model 

We will compute the anomalous magnetic moment of the muon in a general = 1 
globally supersymmetric model containing a U{1) gauge vector superfield, V = (A, A^, D), 
and three chiral superfields, M = {fl,fi,F^^), = (/I^, /i^ F^'-^'), Z = {z.ip^^F^), with 
the following charge assignments: Q{M) = —1, Q{M'^) = +1, Q{Z) = 0. Despite its 
simple field content, the model should reproduce all the relevant aspects of the realistic 
case. In particular, the U{1) will be associated with the exact gauge symmetry of su- 
persymmetric QED, the M and M'^ multiplets will contain the degrees of freedom of the 
left-handed muon and anti-muon, respectively, and the Z multiplet will contain the gold- 
stino. The most general effective Lagrangian with the above field content is determined, 
up to higher-derivative terms, by a superpotential w{Z, M, M'^), a gauge kinetic function 
f{Z,M,M'') and a Kahler potential K{Z,'Z, M,M, M",!^). We can parametrize such 
functions as follows: 

w = FZ + -Z^ + niMM^ + pZMM" + . . . , (3) 
6 

1 / 2Z MM" \ 



|4 



\Z\ 
4A2 



|Z|2|M|2 |Z|2|M^|2 (Z^MM'^ -f h.c. 

"'^^ 2A^ + ■ ■ • • 
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In the above expressions, all the parameters have been taken to be real for simplicity, and 
the dots stand for terms that either do not play any role in the calculation of or can be 
eliminated by (non-linear) analytic field redefinitions^. For example, a possible cubic term 
in K of the form {ISkZMM^ / K + h.c.) can be eliminated by the analytic field redefinition 
Z [Z — PkMM'^/A), accompanied by the parameter redefinitions (m — /SkF/A) —>■ m 
and (7j — ^tjiSk) —>■ 7/- Other field redefinitions have been used to eliminate other low- 
order terms and to shift (z) to zero. Taking all this into account, the above expressions 
of w, f and K are the most general ones compatible with a classical vacuum with broken 
supersymmetry and unbroken CP and U{1). Indeed, the parametrization is still slightly 
redundant, since one of the dimensionless parameters, e.g. rj, can be re-absorbed in the 
definition of the mass scale A characterizing the non-renormalizable operators, but we 
stick to it to have a clearer physical interpretation of our results. 

By standard techniques, it is easy to check that there is a local minimum of the classical 
potential where supersymmetry is spontaneously broken, with vacuum energy {V) = F^, 
and the gauge symmetry remains unbroken, with the gauge coupling constant given by 
((Re/)~^) = e^. Notice that the Kahler metric is canonical at the minimum, so that the 
component fields of the chiral superfields are automatically normalized. In the fermion 
sector of the model, ipz is the massless goldstino, whereas the muon has a Dirac mass 

= m, and the photino has a Majorana mass Mx = rjF/A. In the scalar sector of the 
model, the real and imaginary parts of the complex sgoldstino, z = (S + iP) / -\/2, have 
squared masses: 

ml = az— + aF, w?p = a^— - aF , (6) 
and the smuons have the following squared mass matrix: 

Mi=(<;< ,V (7) 



^2 ^ ~2 ^ r_2 



where 

We denote the smuon mass eigenstates by fli = cos 9fi + sin6'/2'^, /i2 = — sin 9 fl + cosOjl'^, 
and the corresponding mass eigenvalues by m\ and m|, respectively. We also recall that 

sin2^ = ^^ ^. 9 

mf — 1712 

In the model defined by eqs. the different classes of one- loop Feynman diagrams 

that may contribute to are displayed in fig. ||. The Feynman rules needed to compute 

^The above parametrization corresponds essentially to choosing normal coordinates |]l3|] , and simplifies 
our diagrammatic computations. Although reparametrization covariance is lost in the intermediate steps, 
general arguments guarantee that the final physical results are the same. We have explicitly checked this, 
starting from a fully general parametrization. 
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all these diagrams are summarized in the following Lagrangian 
L = -^F^,F'^'' + i\a'^d^\-^[\\ + \\) + hd'^S){d^S)-^mlS^ 



1 



1 



2 I ~ 1 2 

- m2|/i2| 



72 F 



2^2 P 



4v^ F 



/XI/ ^ per 



— — I (m^ — m^) cos6'/ii — — m^) sin6'/i2 /iV'z + h.c.j 
-— I (m^ - m^) sinOfli + (mg - m^) cos9ji2 l^^^z + h.c.j 



+ev^ (cos 6'yUi — sin 6'/i2)/iA — (sin 6'/ii + cos 6jj,2)l^^^ + h.c. 
-7/^^^ {tflli''(T^''\F^u + iil^ixa^"'\F^,, + h.c.) + . . . , 



(10) 



where we used two-component spinor notation in the conventions of ref. [|T3| (with canon- 
ically normalized photon and photino fields, and the field redefinition A iX), and the 
dots stand for the gauge-fixing term and other terms that are not relevant for the following 
calculations. In eq. (10), e°^^^ = —60123 = — 1, and the covariant derivatives are defined 



according to the following sign convention: = + ieQA^. 

To conclude the description of the model, we observe that almost all the parameters 



controlling the interactions of eq. (10) are related with the supersymmetry-breaking scale, 
a/F, and the particle spectrum (masses and mixing angles): the only two additional 
parameters are 7/ and 7^. 



3. General results 

We now move to the computation of the different diagrams depicted in fig. 1 and to 
the presentation of the results, without any assumption on the model parameters. 

The diagrams in figs, la and lb are nothing else than the muon-muon-photon and 
smuon-smuon-photino diagrams of supersymmetric QED, whose contributions are well- 
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known. Fig. la gives the one- loop QED result [p!5[] : 



2n 



(11) 



Its supersymmetric counterpart, fig. lb, gives 



a 2 
— m„ 

271 ^ 



dx 



X 



X] 



a dm? Mx 



m\x + M|(l — x) — m^x(l — x) 
x{l — x) 



+ (1 



7r mf 



ml 



^^mfx + M|(l 



X] 



m^x(l 



'1 



(12) 



Most of the remaining diagrams of fig. |T| involve goldstino or sgoldstino exchanges. 
Some of them are finite, others lead to potential divergences and require a regulariza- 
tion. We will use two alternative regularizations, appropriate for supersymmetry: either 
dimensional reduction in D-dimensional momentum space, or a naive cutoff Kuv in 4- 
dimensional momentum space. To keep track of the regularization dependence, we will 
use the abbreviations: 



A 



uv 



jz^ — 1E + log47r (dim. reduction 
log — ^ — 1 (mom. cutoff ) 



A 



1 (dim. reduction 
(mom. cutoff ) 



(13) 
(14) 



where je is the Euler-Mascheroni constant and fi is the renormalization scale. 

The diagrams of class (c), containing smuon-smuon-goldstino loops, give a finite result: 



m. 



167r2F2 



m2)2x(l — x) 



1 [mf 

dx 

mf — (1 — x) 



+ a 



(15) 



The diagrams of class (d), containing photino-smuon-goldstino loops, are superficially 
divergent. However, the overall result is again finite: 



87r2F2 



A + 



2ml 



mlMl 



■1 rl—x 

„ r, .dx dy 
mf — m2 Jo Jo 



ml)x 



[mf — m2 )x' 



mt 



1 rl-x I 

„^ , dx dy { — 5- „ 

Jo Jo \mfx + Mxy — m^x(l — x 



mfx + M^y — m^^xil — x) 
+ (1-2)) , 




(16) 



where A was defined in eq. (p!^. Notice that the first term within brackets is regularization 
dependent: we will see shortly that it is cancelled by a similar term from another diagram. 
The contribution of the photino-smuon diagrams of class (e) is in general divergent: 

/•I / m?a; + M?(l — x) — mf,x(l — x) , ^ 
-7/T7-^ [2A^v - (log + (1 - 2)^ 

1 /■! m\x ^ M\{\ — x) — m\x{\ — x) 

' dx X log ■ 



-If 



167r2A2 
mf,6m'^ 



87r2A2 m\ — m^ Jo 



mix + M?(l 



x) — m2x(l — x) 



(17) 



where /S.uv was defined in eq. (|I3|). If 7/ = 0, tlien 5a(f) = 0. 

Tlie diagrams of class (f), containing muon-muon-sgoldstino loops, give a finite result: 



167r2 



+ 



5m? 


F 


F 




5inr? 




F 





x^{2 - x) 



10 m|(l — x) + m^x^ 

dx 

mp(l — x) + irfix'^ 



The overall contribution of the diagrams of class (g), containing the photon-muon- 
sgoldstino loops, is in general divergent: 



5af 



IK 



m^Mx 



47r2A2 



A 



uv 



1 
2 



1 

dx 



l-x 



dy [ log , ^ +(g ^ P)4 



2 2 
m,,x 



mi,5m?M\ 



8n^F^ 



A - 2 /dx 



, M mly + mlx'^ 

dy I log -zi 



2 2 
m,,x 



rripy + m^x^ mpy + m^x^ 



(19) 



where Af/y and A were defined in eqs. ([T3| ) and (0), respectively. If = 0, the result 
is finite. Notice also that the regularization-dependent contribution in the second line of 
eq. (|19|) cancels exactly with the one of eq. (|I6|). 

As a consistency check of our results, we have verified that vanishes in the limit of 
exact supersymmetry, as expected on general grounds We recall that, in the unbroken 
phase of supersymmetric QED, the overall is zero, as we can easily check from eq. (12) 
by putting ml 



m% 



and M\ 



a 
'2^ 



(20) 



Similar cancellations take place among the other diagrams. Indeed, if we re-express the 
spectrum parameters in terms of the original ones and take the appropriate F —>■ limit, 
we find: 



6af 



Tjpm 



(A + 1) 



0. 



(21) 



Incidentally, this confirms that the cancellation of the terms ±m^6m'^MxA/[8TT'^F'^) in the 
general results above is just a consequence of supersymmetry. 



4. Discussion 

We now discuss the interpretation of our results. We start from the observation that, 
in an effective theory with spontaneously broken global supersymmetry and arbitrary 
defining functions {w,K,f), the one-loop contribution to is in general divergent [see 
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eqs. (^) and (|T9|)]. This is in agreement with the results of ref. |]T0[, obtained within the 
supergravity formahsm. It is also consistent with the fact that we could have supplemented 
our standard (two-derivative) effective supersymmetric Lagrangian with higher derivative 
operators contributing directly to a^. An example of such a 'counterterm' is 0] 



c 

A3 



ZWM M" + h.c. 



(22) 



where is the supersymmetric field strength, the covariant derivative of supersym- 
metry and c an arbitrary coefficient. The expansion of (^) in component fields contains 
the operator 

2?p 7 7 

-— FV^V^M^ + h.c, (23) 

which in turn becomes a magnetic moment operator^ if supersymmetry is spontaneously 
broken by a non- vanishing VEV for the auxiliary field . 

It is useful to reinterpret the divergent contributions to found in eqs. (^) and 
([T9| ) using a similar language, separating the perturbative generation of supersymmetric 
operators from the question whether supersymmetry is broken or not. Thus we could say 
that our effective supersymmetric Lagrangian generates an operator of the form (|23|) at the 
one-loop level, with a logarithmically divergent coefficient c = 77(7/ — 27/f) A[/y/(327r^). We 
can go further and try to characterize such divergences in a more general way, manifestly 
invariant under analytic field redefinitions. To this purpose, we observe that c itself should 
be viewed as the VEV of a field dependent, reparametrization covariant expression, i.e. 
the operator should be written as0 

where Di {i = Z, M, M'^) is the scalar field reparametrization covariant derivative and 
R\^^ the curvature tensor of the Kahler manifold: 

D,D,f = - fiiK-'Y^K^,, , R\^^ = [K-'Y^iKrmnj - K^p{K-y~m-,,j) . (25) 

For example, with our coordinate choice of eqs. (0)-(0), only the 'jfMM'^ term in / and the 
'JkZ'^ MM'^ term in K contribute to {DMDM'^f) and (-R^^^^^^^), respectively. We recall 
that other terms that could contribute to {OmDmc/) and {R^]^,jzm<=) ^^re removed by 
suitable field redefinitions. Our general expression (|2^ ) can be also used to reinterpret 
the divergent results of [jl^], where no MM'^ term was included in /, but a ZMM'^ term 



was present in K. We can view this as an alternative coordinate choice, which leads to 

^ For convenience, in this section F^n denotes the original (not normahzed) gauge field. Also, we use 
the definition = Z\qq, so that (F^) = — F in our model. 

^The scalar field dependence of the fermionic terms that we are discussing is similar to that of certain 
bosonic terms of ref. where the logarithmically divergent bosonic contributions to the one- loop 

effective action of supergravity models were computed. It is plausible that our terms are supersymmetric 
partners of some of those of ref. pq|, although we were unable to establish such a connection in full detail. 
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a non- vanishing (DMDM'^f) through the connection term [more specifically, from a type- 
(d) rather than a type-(e) diagram]. We should add that the supergravity calculation 
of ||T0[ found several divergent contributions in individual diagrams, most of which were 
cancelling in the final result and two of which were not. In our case, the use of a globally 
super symmetric Lagrangian (in a convenient field parametrization) avoids the proliferation 
of divergences and reproduces directly the two 'genuine' ones. 

Our covariant formula (|^) clearly shows the geometrical meaning of the two kind of 
divergences. We could now ask whether they should be considered as independent or could 
possibly originate from a single object. In fact, taking into account that {fu) = (/a/=) = 0, 
we observe that the VEV in eq. (H) is identical to the VEV of —AD-^DMDM'=^ogRef. 
Thus we are led to conjecture that a (two-derivative) supersymmetric Lagrangian like 
ours, with generic / and K functions, generates radiatively the following logarithmically 
divergent operator: 

{D^DmDm^ log Re/) F^fia'^y^F,, + h.c. . (26) 

The corresponding (higher-derivative) superfield operator should be an appropriate gen- 
eralization of the operator (|2^) . 

We have discussed so far the structure of the supersymmetric operators that give 
divergent contributions to once supersymmetry is broken. Since the existence of such 
operators limits somewhat our predictive power, we could look for models where the 



corresponding coefficients vanish, or are at least suppressed. The simplest case [|, ^, |T0[ 
is to have (fz) = 0, implying Mx = 0. More interestingly, the whole combination in 
eq. ( P^ [or the corresponding compact expression in eq. (PB])] might vanish if / and K 
were related, either functionally or just on the vacuum (the latter case corresponds to the 
equality 7/ = if we use our original parametrization). This approach may lead to 
models with an underlying extended supersymmetry. 

The question could also be addressed in terms of ordinary symmetry arguments. We 
may ask whether there is a symmetry that allows for non-vanishing values of and 
but forbids the divergent one-loop contributions. A possible candidate is the R-symmetry 
defined by the following charge assignments: R{Z) = +2, R{MM'^) = +2 (the separate 
values of the R-charges of M and M'^ are not important here). Such an assignment, 
however, would put to zero not only the coefficients 7/ and 7^^ of the divergent operators 
(and possible coefficients related to them by analytic field redefinitions compatible with 
the R-symmetry), but also the mixing terms in the smuon and sgoldstino mass matrices 
and, most importantly, the photino mass Mx. With the latter constraint, the construction 
of a realistic model with the full Standard Model gauge group appears very unlikely. We 
were unable to find alternative symmetries that forbid the divergent one-loop contributions 
to in the presence of a reahstic spectrum. 

A milder requirement may be to ask for a symmetry that forces the divergent contri- 
butions, eqs. ( P!7D and (|T^, to be at least proportional to m^. If so, reasonable choices 
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of the ultraviolet cutoff would produce contributions to that are at most of the same 
order of magnitude as the finite ones. An obvious candidate for such a symmetry is a 
chiral U{l)s, with the charge assignments: S{Z) = 0, S{MM'^) ^ 0. Such a symmetry 
would be explicitly broken by a small parameter e to allow for the muon mass term. As 
a consequence, also 5m^, 7/ and would be suppressed by the same small parameter. 
In particular, we would obtain ~ m^A, with A of order the other supersymmetry- 
breaking masses. Such a hierarchy would be stable against radiative corrections, as can 
be easily checked by looking at the one- loop diagrams contributing to and to 5rn?. 

In order to discuss in more detail the relative size of the 0{m?^ contributions to a^, we 
should supplement the chiral power counting with additional information. For example, 
we could assume^ that F/ h? ~ e' is small and consider an expansion in both e and e' . 
Even so, at least two possibilities could be considered for the combined power counting: (i) 
•jK ^ If ^ e, ~ eA ; (ii) ■yj^ jjr e, ~ ee'A . In case (i), the contributions to 
proportional to jk and 7/ [or induced by the counterterm of eq. (p2D] would be suppressed 
by a factor e' with respect to the remaining ones, and could be considered as subleading. 
Then we would obtain an unambiguous result at the leading order. In case (ii), on the 
other hand, such contributions would be formally of the same order of the other ones0. 
Moreover, as we shall discuss in more detail below, phenomenological considerations will 
lead us to consider, among others, values of F/A^ of order one. 

We can summarize the above discussion as follows. The effective supersymmetric 
low-energy theory can contain higher-derivative terms that contribute to and do not 
depend (only) on the spectrum. Such terms could either be present as counterterms 
since the beginning, or be generated radiatively within the low-energy theory itself, or 
both. Symmetry considerations could help in keeping the size of such contributions under 
control, as in the example of a chiral symmetry, possibly combined with an expansion in 
F/ h?. In most phenomenologically relevant situations, however, it seems that only a better 
knowledge of the underlying microscopic theory could remove the residual ambiguities. 
Taking into account all this, in the following we will use the finite contributions to only 
to derive some 'naturalness' constraints on the model parameters, barring the possibility 
of accidental (or miraculous) cancellations. 



^We thank the referee for raising this point. 

^Notice that the doubly suppressed muon mass of case (ii) would naturally arise from a Kahler potential 
term of the form PkZMM'^ / A., with fix ^ £■ More generally, case (ii) can be justified by an R-symmetry 
under which Z, M, M'^ are neutral, explicitly broken by the parameter e': then all superpotential parame- 
ters would be suppressed by e'. Incidentally, this would also imply the same scaling law for all sgoldstino 
mass terms, (m| ± mp) ^ F^/A^. 
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5. Phenomenology 



Keeping in mind the discussion of the previous section, we now focus on the class of 
models in which 

{DMDM^f) = , {Rl.^.J = 0, (27) 



so that the total one-loop contribution to is finite and determined only by the param- 
eters controlling the spectrumQ. We specialize the general formulae of sect. 3 to some 
interesting limits, and we discuss the resulting phenomenological constraints. 

We begin with the SQED contribution. For realistic smuon and photino masses, so 
that the muon mass becomes negligible in the comparison, eq. (|T2|) becomes 



a 2 
m,, 

271 ^ 



a 5m? M\ 
n nil ~ "^2 



Qm\Ml + Sm^M^ + 2Ml + Gm^M^ \og{m\/Ml) 



Mi 



6(mf - MlY 
2mlMl \og{ml/Ml) 



+ (1-2) 



2(mf 



Mlf 



(28) 



Next, we consider the contributions from the smuon-photino-goldstino sector, for 7/ = 
0. Eq. (p!?!) gives no contribution, whereas eqs. ( plSf ) and ([I^) simplify considerably in the 
phenomenologically relevant limit <^ mi, 1712, Mx: 



m'j^{m1 + m"^) 
967r2F2 



m,.6m'^Mx 
^ -A + 



m. 



167r2F2 



mjMl 



Ml 



+ (1 



(29) 



(30) 



m{- Mi \ 

The result in eq. ( P^j ) is in full agreement with ref. and disagrees in sign with the 
supergravity computations of refs. |TD[. The result in eq. ( PI1| ) can be compared with a 
similar expression of ref. |TD[ , where the terms within brackets appear with opposite sign. 
Also, we recall that the first term of eq. (^) is regularization-dependent and disappears 
when all contributions are summed. 

We now discuss the contributions from the muon-photon-sgoldstino sector, eqs. (|18[) 
and (|19D, for 7^^ = 0. Before taking any kinematical limit, we would like to emphasize 
that such contributions depend only on the spectrum parameters and the supersymmetry 
breaking scale, with no more model-dependence than the contributions from the smuon- 
photino-goldstino sector. This point was apparently overlooked in ref. [n)|. However, we 
have checked that the results for the muon-photon-sgoldstino sector given there agree with 



^The finiteness conditions ( [27| ) generalize those of ref. |T^. In our original parametrization, they 
amount to set 7/ = = 0. We will not discuss here the alternative intriguing possibility mentioned 
in the previous section, corresponding to 7/ — 7^ 0. We only recall that also in the latter case 
the result would be finite and regularization-independent, as long as we use regularizations compatible 
with supersymmetry. However, it would contain an additional parameter besides those controlling the 
spectrum. 
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ours, once the flat limit is taken and appropriate parameter identifications are made [e.g., 
3m^/2{mf^ — m") = Sm"^]. Furthermore, the general results in eqs. ([T8|) and ([T9|) simplify 
considerably in the limit of heavy or light sgoldstinos: 



H"-+S^H .(31) 

c- , ^ c- o , ^ f log ^ imn^mp.mQ) 

'^^'^ "+ 87r2F2 ^+ 87r2F2 • ^^^^ 

[ 1 (mp, < m^) 

Taking other special limits such as mp(^s) ^ "^/^ "^siP) is equally straightforward. We 
recall that the possibility of superlight sgoldstinos (both of them or just one) has been 



frequently considered in the superlight gravitino literature [17]. Despite its naturalness 



problems []T8l , it may be related with possible dynamical mechanisms for the resolution of 



the strong-CP and/or the cosmological constant problems. 

We conclude this section by confronting the above results with the present and future 
experimental constraints^. In the following, it will be convenient to parametrize the off- 
diagonal element of the smuon mass matrix as Sm? = m^A, as suggested for example by 
an underlying approximate chiral symmetry (see sect. 4). Thus, although A, as defined 
by the previous relation, is a free parameter equivalent to 5m^, it is natural to expect 
that A be not much larger than the other supersymmetry breaking masses. We begin 
by reconsidering the SQED contribution of eq. (|28D . We recall that such a contribution 
should be replaced, in a fully realistic model, by the full MSSM contribution |^: here it 
will be used only as a benchmark for the new contributions associated with the superlight 
gravitino and its superpartners. For illustration purposes, we can further simplify eq. ( pSf ) 
and consider just a representative case. For example, if the smuons and the photino are 
roughly degenerate, with a common mass M5, we obtain 

c (h\ a ml / A\ in (100GeV)2 / A\ 

I" 127rM|V Ms J Ml \ Ms) ' 

As we can see, sizeable contributions can be obtained only if the sparticle masses are 
significantly smaller than 100 GeV or if the supersymmetry-breaking mass parameter A 
is significantly larger than the typical supersymmetry-breaking masse^. 

^For previous analyses based on the one-loop goldstino contributions of ref. |T^, see refs. jlll and 
[ p^ . The latter paper studies also the influence of superlight sgoldstinos, but neglects their effects at the 
one-loop level and only partially includes them at the two-loop level. 

^We recall that, in the framework of the MSSM, the role of our effective A parameter would be played 
by the MSSM parameters A^^ and /itan/3. Thus, for example, enhancement effects could originate from 
large values of tan /3 |7j . 
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To discuss the phenomenological impact of the remaining contributions, it will be useful 
to parametrize all of them in a uniform fashion 



(34) 



where the real parameter has the dimension of a mass squared and can be positive or 
negative. The above form separates the dependence of 5a ^ on the spectrum (through Mj.) 
from that on the supersymmetry-breaking scale a/F. Incidentally, we recall that the latter 
scale is related to the gravitino mass r?7,3/2, the reduced Planck mass Mp and Newton's 
constant Gjq as follows: = 37713^2^!- = 37713/2/(8716' at). The present experimental limit 
and the future sensitivity on are then shown in fig. ^ the two solid lines correspond to 
the present upper and lower bounds on 5a^, see eq. (|]), whereas the dashed line corresponds 
to \5a^\ = 4 X 10^^°, and gives a measure of the expected future sensitivity. 



1000 



800 



600 



> 
a 




1000 



800 



> 
a 



400 — 



200 — 



600 



400 



200 



200 300 
VF (GeV) 



500 




I I I I I I I I 



200 300 
VF (GeV) 



400 500 



Figure 2: Contours of 5a^ = rrr^M^ / [IQix F ) in the plane (v-F, |Mj;|).- the regions above 
the solid lines correspond to 5a^ < —70 (> +220) x 10~^°, and the dashed lines correspond 
to 5a^ = +4 X 10-1°. 



The information in fig. 2 should be combined with the explicit expression of in 
terms of the spectrum, which can be easily read off our results. In particular, we can 
put together our results of eqs. (p9|)-(|3^) and express as the sum of 'goldstino' and 
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'sgoldstino' contributions, 

Ml = Ml^a) + Ml^sG) , (35) 

where 

M\g) = < 

and 

r 2AMx log ^ + 0{A^ml/m\p) (m^ < mp, ms) 
M.'(5G)= ' • (37) 

{2AMx + A^ {ms.mp <^m^,) 

Notice that we have aheady removed from each contribution the regularization-dependent 
terms that disappear in the overall result. 

The size of the goldstino contribution to depends mainly on the smuon masses, 
and can be significant only for heavy smuons and a very low supersymmetry-breaking 
scale. The corresponding sign changes from positive to negative for increasing smuon-to- 
photino mass ratio. If the smuons and the photino have roughly the same mass M5, we 
obtain M^^^^ ^ 2M|/3. In interpreting fig. |^, however, we should keep in mind that the 
situations where |Ma.(G')| ^ ^/F correspond to a strong-coupling regime, and are limited 
by perturbative unitarity. As discussed for example in |T^, [T^ , a reasonable requirement 
is |M,(G)| < (2-3)v^. 

The size of the sgoldstino contribution to depends on the photino mass, the sgold- 
stino masses and the parameter A. The latter parameter, which controls the off-diagonal 
element in the smuon mass matrix, plays a crucial role. In particular, the sgoldstino 
contribution is zero for vanishing A. Several possibilities can arise for non-vanishing A. 
We comment here on some representative cases, using the symbol Ms to denote a typical 
smuon or photino mass. Intermediate cases not discussed below can be easily worked out 
from our general results. 

i) Heavy sgoldstinos. By 'heavy' here we mean the case msyfrip = 0{Ms), which is 

favoured by naturalness considerations [0. Then from eq. ( ^7\ ) we have M^^^g^^ ^ 
2^Ma log (mp/m|), with either sign being possible. Notice that the sgoldstino con- 
tribution becomes negligible if S and P are approximately degenerate. Otherwise, 
we should also specify the size of A. In the natural case of A = 0{Ms), M^^^q is 
0{Mg) as well, thus the sgoldstino contribution to is C(m^M|/167r^F^), like the 
goldstino one. However, if (for some reason) A 3> Ms, the sgoldstino contribution 
gets enhanced by a factor 0{A/Ms) compared to the goldstino one [for a similar 
situation in SQED, see eq. (PB|)]. 

ii) Superlight sgoldstinos. By 'superlight' we mean here the case ms,mp -C (however, 

the qualitative picture remains the same whenever at least one of the sgoldstino 
masses is < m^). Then from eq. ( ^7]) we have Ml,gQ^ = 2AMx + . Again, if 



Ml 



Ml 



6 m\-Ml 



log- 



m\ - Ml ^ Ml 



2) 



(36) 
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A = 0{Ms), then M^^g^^^ is 0(M|) (either sign being possible), and the sgoldstino 
contribution is comparable with the goldstino one. If instead A ^ Ms, then we have 
M^f^gQ-^ ~ > Ml^fjy and the sgoldstino contribution to (= m^A^/ievr^F^) 
dominates over the goldstino one. Thus the simultaneous occurrence of superlight 
sgoldstinos, large A parameter and low supersymmetry-breaking scale gives a po- 
tentially large contribution to a^, comparable with (or even exceeding) the present 
experimental bounds. However, it should be taken into account that such a situ- 
ation is an extreme one, as we have already remarked. Moreover, other potential 
enhancement effects due to superlight sgoldstinos have been studied in the context 
of cosmology, astrophysics or collider physics jl^, and the inferred lower bounds on 
\fF are already well above the TeV range, weakening the impact of the constraint. 

In summary, since the present lower limits from accelerator searches should be of the 
order of 100 GeV for the smuon and photino masses[^ and of the order of 200 GeV for \/F 



20| , we can see that at present provides a non-negligible but mild indirect constraint. 



Such a constraint, however, will become much more stringent after the completion of the 
E821 experiment. Should a discrepancy emerge between the future E821 result and the 
SM prediction, models with a superlight gravitino may provide a viable explanation. 



6. Final comments 

We conclude this work with some comments on the predictive power of our effective- 
theory approach, to clarify its advantages with respect to non-linear realizations of super- 
symmetry, as well as its residual limitations. 

Since is a low-energy observable involving only SM particles, we may naively think 
that its calculation could be addressed in the context of a 'more effective' theory, where 
heavy particles have been integrated out and supersymmetry is non-linearly realized. In 
the case of heavy sgoldstinos, a first step in this direction would be to integrate out only the 
sgoldstino fields, ending up with a non-linear goldstino superfield I^Tj . A more radical step 



would be to integrate out not only the sgoldstinos but also the smuons and the photino, 
ending up with a fully non-linear realization of supersymmetry on the muon, photon 
and goldstino fields. It is immediate to realize, however, that non-linear realizations are 
not a sufficiently predictive approach to address the computation of a^. This should be 
contrasted with other processes involving external goldstinos [|1^, [2^] , where the low-energy 



theorems of spontaneously broken supersymmetry can be put to work. The simplest way 
to see this point is to look at the dependence of our general finite result ( |5^ ) on the 
supersymmetry-breaking masses, rewriting it as: 



""^^A dedicated analysis for models with a superlight gravitino is not available yet, and would require the 
generalization of the model to include the full SM gauge group. 
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where Ms is a typical supersymmetry-breaking mass that can be taken as a reference 
scale, the m? are the individual supersymmetry-breaking mass parameters in the sgold- 
stino, smuon and photino sectors, and the dimensionless function G encodes detailed 
information on the supersymmetry-breaking spectrum. From the explicit form of G [see, 
e.g., eqs. ([55|)-(P7D], it is evident that there are no decoupling properties for particles with 
large supersymmetry-breaking masses. This result is expected: decoupling should not oc- 
cur for heavy particles whose masses break a symmetry that would set the result to zero. 
In the case of a non-linear realization, we could at most obtain the general scaling law of 
eq. (p8D, but the function G would be just replaced by a free parameter (if we replace the 
linear superfield Z by a non-linear goldstino superfield, also the form and the dimension of 
the lowest-dimensional counterterm for will change). Starting from an effective theory 
where some of the superpartners are not present would in general reduce our predictive 
power, generating new divergences without a prescription on how to regulate them. In 
this respect, removing only the sgoldstinos, while keeping the smuons and the photino, 
would not lead to any substantial increase in predictive power: sgoldstinos do not have 
any special decoupling properties with respect to other particles that may acquire large 
supersymmetry-breaking masses. 

The second and last comment is a warning on the correct interpretation of the phe- 
nomenological implications of our results. As already discussed in sect. 4, the finite result 
of our general calculation, obtained by setting ■Jk = If = 0, is protected against large cor- 
rections only when a/F <^ A and there is a symmetry guaranteeing that 7^ ~ 7^ ~ m^/A. 
In this case, the only realistic possibility of generating observable contributions to is 
to have A much larger than the other supersymmetry-breaking masses, to enhance the 
sgoldstino contributions of eq. (|37|) . This is certainly possible, and is the exact counter- 
part of having large values for and/or fi tan/5 in the MSSM. However, the value of 
A cannot be pushed too far, since A breaks the same chiral symmetry that protects the 
smallness of m^: too large values for A would make unnaturally small. Apart from the 
above possibility, observable contributions to 6a ^ can be obtained only when 7^^ and 7/ are 
unsuppressed, in which case we have no predictive power at all, or when ^/F ~ A. In the 
latter case, we can still make use of our finite result, but we should keep in mind that we 
are neglecting other contributions of comparable size. If no deviations from the SM pre- 
dictions are found, we can still use the results to derive some 'naturalness' constraints on 
the model parameters, barring the possibility of accidental (or miraculous) cancellations. 
Should a discrepancy emerge between the experimental data and the SM prediction, how- 
ever, a detailed quantitative analysis would require a better knowledge of the underlying 
microscopic theory, to remove the ambiguities associated with the operators controlled by 
jK and 7/ and with possible unsuppressed higher-derivative supersymmetric operators. 
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